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Abstract 

We suggest a simple analytical description of the S-wave isoscalar tttt amplitude, which 
corresponds to a joint dressing of the bare resonance and background contributions. The 
amplitude describes well the experimental data on the 5® phase shift in the energy region 
below 900 MeV and has two poles in the Re s > half-plane. Besides the well-known pole 
of <7(600)-meson with Re s ~ there exists a more distant pole with Re s ~ 0.6 GeV 2 . 
Our analysis indicates for the dynamical origin of the cr(600) pole, while the second pole 
should be associated with lowest qq state. 

~cn~ 
o . 

^ : 1 Introduction 

o 

The properties of the lightest scalar meson <r(600) are very important for interpretation of a 
scalar family and details of the chiral symmetry breaking. Appearance of a new experimental 
information and the theory development in the low energy region generated an extensive dis- 
cussion on this issue (see [TJ El El HI E] and references therein). As for existence of <r(600), now 
it is a commonly accepted fact and this resonance is included back into Particle Data Group's 
tables. 

There is a long story concerning the resonance interpretation of the S-wave amplitude tvtt — > 
7T7T with isospin 1=0. One of the key moments of this story was a realization (see e.g. [HI IB!) 
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that apart of the resonance a (600) term there is an essential background contribution in this 
energy region. However there is no evident recipe to divide the tttt amplitude into the resonance 
and background terms. The simplest and widely used method is " adding in phase shift" of the 
o ■ resonance an background contributions (IA method in terminology of 0IE]) : 

6° = 6 R + 6 B . 

^ ■ The amplitude tttt — > tttt in this case is: 
& • 

„2i<5° i „2i<5 s i p 2i5 R i 



2ip 2ip 2ip 

+ e 2iSB ■ f Res (s). (2) 



e 2iS« _ ! 



2ip 

The anzats ((IJ may be derived from summation of the loop contributions with some extra 
conditions jHj. To obtain the resonance parameters from the experimental data one needs an 
additional assumption about the form of the background contribution 8 . The best way for a 
broad resonance is to determine its mass and the width from the pole position in the complex 
energy plane. However only the resonance contribution f Res (s) of the entire amplitude (j2J) can 
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be continued into the complex energy plane. Thus it is possibly to study the pole position but 
not the pole residue. 

Other methods to describe the 5® phase shift, different from (pQ), either have the so evident 
defects, or are much more complicated with many free parameters. 

In this paper we suggest a very simple analytical parameterization for tttt amplitude which 
allows us to continue it to the second Riemann sheet. The amplitude contains (in spirit of the 
linear cr-model) two bare objects: the resonance and the background. The main idea is that 
a joint unitarization of two objects should be described correctly by the field theory methods. 
As concerning the form of a background contribution at the tree level, it can be modelled by a 
maximally simple method. 

There are different ways to construct such analytical amplitude. We found the suitable 
one the formalism of the unitary mixing, the obtained amplitude is analytical and unitary 
automatically. Such construction is rather flexible which allows us to investigate some different 
physical situations. 

Note that from other side the bare pole located at s < may be considered as some effective 
cross exchange and its value m 2 , ~ -m^, obtained from a fit, confirms this interpretation. As 
compared with standard N/D method our amplitude with two bare objects automatically has 
a zero, which is necessary to describe the S-wave low energy data. 



2 Formalism of the unitary mixing 

If there exist n bare states with the same quantum numbers then the dressing of their propa- 
gators should account also the mutual transitions between them. The process of joint dressing 
is described in this case by the system of Dyson-Schwinger equations: 

TLij(s) = Trij(s) - U ik (s)J k i(s)7rij(s), i,j = l...n. (3) 

Here tt^ and are bare and dressed propagators respectively, Jij are the self-energy 
contributions. 

Let us consider mixing of two resonances (n=2) with one open intermediate state. In this 
case the solutions of © are: 

n -° 2{s) n - Ms) n - Dl{s) U) 

Here 

D(s) = D l D 2 -{J 12 {s))\ 

Di = m\ — s — J n {s), D 2 = m 2 2 - s - J 22 (s). (5) 

In the case of scalar resonances interacting with pion pair, the loops are of the form: 2 

Jn(s) = gfJ(s), J 22 (s) = g 2 J(s), J 12 (s) = gig 2 J(s). (6) 

2 Note that we ignore the subtraction constants in the loops. As it is shown in Appendix A a subtraction 
polynomial in the loops can be removed by the redefinition of bare parameters. 
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. s - a f°° ds' ... . . Is - 4ml . . 

J(s) = — L p(s ) ' p(s) = v — r- 1 - (7) 

where < a < 4m 2 is the subtraction point, ^ are the coupling constants. 
The 7T7T amplitude: 

/ = glu n ( s ) + g 2 2 u 22 ( s ) + 2< M2 n 12 ( s ) = (8) 

D{s) 

where 

D(s) = (m{- s)(m 2 - s) - J(s)N(s), 

N(s) = g 2 1 (m 2 2 -s)+g 2 2 (m 2 1 -s). (9) 

It is evident that (JBJ) satisfies the elastic unitary condition 

Imf = p\f\ 2 . (10) 

The above equations can be applied not only for the case of two resonances but also for the 
"resonance+background" situation, when one of the bare poles is located at s < 0. Just right 
this situation arises for the S-wave 1=0 tttt — > tttt amplitude. One can see from (j!2)) that our 
amplitude is zero at the point: 

So = (9iml + g 2 2 ml)/(9 2 i+9 2 2), (H) 

which should be s[j ~ m 2 to reproduce the Adler zero. So we have m\ < 3 . 

Let m\ < be the bare zero of function D(s), which stays at the left of real axis after 
dressing. Then it is convenient to subtract the loop at this point: 



D(s) = (m 1 - s)(m 2 - a) - (J(s) - J{m 2 )) g x (m 2 - s) + g 2 (m 1 - s) 



'12) 



Below we shall use the amplitude (|H|) . ()12jl for description of the experimental data. Here 
m 2 , gf, m|, g\ are free parameters. The background contribution at the tree level may be mod- 
elled by the pole or constant. It is sufficient for successful description of the experimental data 
as it is seen below. 



3 Analysis of tttt data in region of < 900 MeV 

In the nearthreshold region we use the new data from decay ^0], which may be seen at 
Fig. ^ in comparison with 1977 year data JT]. We do not take into account the old data jTT] as 
it has no practical effect on the fit. The measured value in decay is the phase shift difference 
Sq — 5{ thus we need an additional information on the P-wave. We use for this purpose the 
approximation of solution of the Roy equations from Ref . ^2] • Fortunately, the 5{ contribution 
is only about 1.5° at the end of the interval due to the P-wave threshold behavior, thus the 
uncertainty in 5{ is negligible. 

3 In spite of < 0, we keep using this notation to stress the presence of two objects in the amplitude. 
Note by the way that our amplitude (JHJ coincides except of notations with the amplitude of Ref. :9; , obtained 
from the low-energy bootstrap equations. 
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Figure 1: Phase shift difference 5q — 5\ from experiment on decay. 



Our main purpose is the a (600) resonance thus we restrict ourselves by the energy region 
m nn < 0.9 GeV. It allows us to use the one-channel approach and not to take into account the 
/o(980) effect. In this region there exist different experiments and different analyses of S-wave 
phase shifts, see recent reviews [131 HI] . 

Below we use only classical partial analyses of Protopopescu et al. ^Hj from 7i + p — > 
7r + /T~A ++ reaction and Estabrooks and Martin one of CERN experiment J2| on n~p — > 7r + 7r _ n 
(let us call their two solutions for 5® phase shift as EM I and EM II). Below we consider the 
mentioned experimental data on the 5{J phase shift and find very similar conclusions. As an 
example we focus in more details on the the EM II solution. 

Fig. El displays the results of joint fitting of data (m nn < 0.4 GeV) and EM II data 
(0.51 < < 0.9 GeV). One can see that our amplitude (|T2"|) describes well these data. 

Best fit parameters are: 



Let us consider the zeros of function D(s) at the second Riemann sheet 4 . The procedure of 
analytical continuation is described in Appendix B. The Fig. El shows the zeros location in the 
complex s plane corresponding to the best fit parameters (fl"3*|) . Let us stress that we find two 
zeros in the Re s > half-plane. The Table ^ represents results of the fit of the different low 
energy data by our amplitude (fT2"|) . All data sets lead to the solutions with two poles: close 
and distant 5 . 

In FigE] we compare phase shifts corresponding to different variants of Table H 

4 The values of the bare parameters have rather limited meaning since they correspond to a given method 
of renormalization. However the character of the pole movement is more meaningful, at least when the loop 
contributions do not dominate in amplitude. 

5 Our amplitude has a property f(s*) — f*(s) thus we have a pair of the complex conjugate poles in the 
complex s plane. For definiteness we say about poles in the Im s > half-plane. 




0.659 ± 0.041 GeV 2 , g\ = 0.435 ± 0.036 GeV 2 , 
-0.230 ± 0.114 GeV 2 , g\ = 0.177 ± 0.067 GeV 2 , 
17.7/21. 
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Figure 2: Results of the fit of the im phase shift by the amplitude (fT2"J) . Solid curve is the best 
fit to the data of K14 decay and EM II its parameters can be seen in the first column of Tabled 
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Figure 3: Poles of the amplitude ttit — > tttt with J=I=0 on the second Riemann sheet at 
parameters (pH?j) . Arrows indicate the direction of the poles movement when the interaction is 
turned off. 



5 



K U +EM II 
E < 0.9 GeV 


K l4 +EM I 
E < 0.9 GeV 


i^/4±Protopopescu 
E < 0.9 GeV 


CGL [12! 

£ < 0.8 GeV 


m\ = 0.659 ± 0.041 
g\ = 0.435 ± 0.036 

~~2 n oon _i_ n i i /i 

wij = —0.230 ± 0.114 
gl = 0.177 ±0.067 
X 2 /DOF = 17.7/21 


mj = 0.586 ± 0.025 
g\ = 0.382 ± 0.020 
m 2 = —0.113 ± 0.053 
g\ = 0.116 ±0.030 
X 2 /DOF = 22.6/21 


mf = 0.794 ±0.114 
gl = 0.598 ±0.151 

~~2 n con _i_ n /inc 

m 2 = —0.580 ± 0.405 
gl = 0.422 ± 0.331 
X 2 /DOF = 5.3/19 


ml = 0.845 
gl = 0.779 

m 2 = —0.573 
gl = 0.548 

X 2 /DOF = 


Poles: 
S! = 0.015 + i 0.192 
s 2 = 0.633 + z 0.630 


Poles: 
si = 0.045 ±i 0.132 
s 2 = 0.632 + i 0.533 


Poles: 
Sl = 0.055 + i 0.339 
s 2 = 0.484 + % 1.020 


Poles: 
si = 0.104 + i 0.250 
s 2 = 0.659 + i 1.620 
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Table 1: Results of the fit of the different sets of experimental data by our amplitude, m. 
Si are in units of GeV 2 . In the last column our parameterization is compared to approximation 
of the phase shift from Ref. [TJ] (Solution 17 of the Roy equations with use of scattering 
lengths from the two-loop chiral perturbation theory calculations.). Our phase shift practically 
coincides with CGL approximation in this energy region. 
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Figure 4: Phase shifts corresponding to dif- 
ferent sets of parameters shown in Table ^ 



Figure 5: Data fit at fixed value of the left 
pole m 2 . 
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We can see that our simple model ()12|) corresponding to joint unitarization of two bare 
objects: one pole at s < and another pole at s > describes successfully the utt phase shift 
5q in the energy region below 900 MeV. We find two poles in the complex s plane: one close to 
the origin with Re s ~ m 2 and the second one with Re s ~ 0.6 GeV 2 . However the behavior of 
the poles when interaction is turned off gf — ^ is rather unexpected (see Fig. |3J): the close pole 
traditionally identified with cx(600) meson moves to the negative s region. While the second 
pole S2 (most of previous analyses did not observe it) tends to the real axis above the threshold. 

As an alternative we can investigate the case when the background has not bare pole. It 
corresponds to the joint dressing in the system "cx-pole + constant". For this purpose it is 
sufficient to put the value m 2 , negative and large in our amplitude (fH2|) . 

In Fig. El there are shown the results of data fit with different m\ values. One can see that 
the experimental data prefer rather close left pole |m|| < 0.6 GeV 2 . 

Fig. IH1 illustrates the pole positions in the complex plane at value m 2 . = — 1 GeV 2 and 
their behavior when interaction is turned off. We observe that the behavior of poles has been 
changed as compared with m 2 . > —0.5 GeV 2 case. 



Re s, GeV 2 Re s, GeV 2 




Figure 6: a) Poles positions in the complex s plane at m 2 = —1 GeV 2 (fit of the data 
K^+EM II) and their movement at g>j — > 0. 

b) Illustration of the poles movement. As compared with a) we slightly reduced the coupling 
constants gf — > 0.75gf with fixed other parameters. 



4 Discussion 

We found that the nir phase shift 5q is well described by a simple analytical amplitude (|12|) 
in the energy region from the threshold up to 900 MeV. Our amplitude corresponds to a joint 
dressing of two bare objects: resonance and background contributions. Background can be 
modelled either by a pole with Re s < or by constant. As a next step one could investigate 
the more complicated background model: left pole + constant (just as in the linear a-model). 
However since our simple amplitude (JT2jl provides a good description of the experimental data 
we suppose that inclusion of new degrees of freedom has no meaning. 
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After the fit of the experimental data we found the presence of two complex poles at the 
second Riemann sheet: one close to the origin with Re s\ ~ m 2 and the second one with 
Re S2 ~ 0.5 — 0.6 GeV 2 . The close pole was seen in most of the previous analyses of ttti 
scattering (its position is defined mainly by Adler zero) and it was associated with the lightest 
scalar meson a(600). Note that we approximated the background term at the tree level by 
some pole, physically it corresponds to the cross-exchange by p- or cr-meson. The existing 
experimental data certainly prefer the pole form of background as compared with constant. 

As for behavior of the poles in the limit of g^ — > we observe that only the distant pole 
goes to the real axis at positive s (see Fig. EJ). This fact holds true for all variants indicated 
in Table E More detailed investigation shows that such a behavior changes with m| value as 
it is schematically illustrated in Fig. The experimental data on tctc scattering with energy 
below 900 MeV prefer the variant a) while the variant b) with m\ ~ — 1 GeV 2 also can not 
be excluded (see Fig. EJ). For example, the found \ 2 values are x 2 /DOF = 17.7/22 for a) and 
X 2 /DOF = 22.6/22 for b) in one of variants of fit. 

Re s, GeV z Re s, GeV z Re s, GeV z 




Figure 7: Poles positions in the complex s plane and their movement at gi — > is managed by 
m 2 value. The variant a) is preferable by quality of the data description but b) also has an 
acceptable value of \ 2 - 

In view of discussion (TBI EE EI3 EU Ej whether the a (600) the intrinsic state or it is 
dynamically generated, our results should be interpreted as an indication for a dynamical 
nature of the <r(600). In this case the second pole should be associated with intrinsic qq state 
having regard to above remarks. 

We suppose that the most interesting question is the meaning of the second pole S2- It was 
seen only in a few previous analysis, e.g. in [2], where it was considered as an artefact since it 
was located out of the considered energy region. In our analysis with account of the much more 
exact data from decay, this pole has moved to lower value Re s 2 ~ 0.6 GeV 2 . As for its 
imaginary part, it is rather uncertain (see Table 1) and may be abnormally large for resonance 
state. We suppose that further fate of this pole may be solved by an analysis in the extended 
energy region. 

In any case it is clear that in fact we have the joint complex " a (600) + Back ground+f o(980)" , 
which should be studied jointly and by the adequate methods. 
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A Reparameterization of hadron amplitude 



Let us consider the unitary mixing of two bare poles with presence of one intermediate state. 
We are interested in a number of independent parameters in tttt — > titt amplitude. Let us write 
it in a matrix form: 

£)(:)-^ n -» (14) 

Here II is the symmetrical matrix of propagator. 

Let us start from the most general case when all loops have a subtraction polynomial of a 
first degree 6 . 

Jij = 9i9j(Pij(s) + J(s)) = g l g j (a ij + bijs + J(s)), (15) 

where 

JW = iz£ r y ,/ZEM. (16) 



7T J4m2 (s' — 0)(s' — s) V s' 

Here a is subtraction point < a < 4m^, for analytical continuation it is not convenient to 
subtract integral at zero. There are ten parameters: bare masses m^mj, coupling constants 
gi, g 2 and 6 subtraction parameters in the loops. 

We can perform a transformation of propagators and coupling constants, which does not 
change the amplitude: 

/ = g T -S- 1 S-U.S- 1 S-g = g' T -U'.g\ 

n' = SUS~\ g = Sg. (17) 
Let us make few transformations consequently: 

1. Firstly by the orthogonal transformation 

S=( COS ° a[n l) (18) 
y — sin v cos v J 

we delete the linear on s term in the non-diagonal loop: b' 12 = 0. 

2. Then by the scale transformation 

we make the coefficient at s in Iln, U22 by unity. After it any orthogonal transformation 
can not generate again the linear on s term in the non-diagonal loop. 

3. We use one more orthogonal transformation to delete a subtraction constant in the non- 
diagonal loop. 

4. Finally, we can redefine the masses, absorbing the subtraction constants in the diagonal 
loops. 

As a result we came to parametrization (fl|)-(|7j) which contains four parameters: masses and 
coupling constants. 

6 Higher degree of polynomials leads to dominating of loops contributions at large s. It leads to the changing 
of the problem's index and to changing of number of poles as compared with the non-interactive case. 
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B Analytical continuation of loop 



Let us consider the two-sheet analytical function: 

V2 e i Vl /2 



Ami 



Ami 



zA ^2/2 



The cuts are chosen from —00 to zero and from Am\ to +00. 

Let us write down the Coshi theorem on the first Riemann sheet 



F«(s)-F«(a) = J (1) +L(s). 



< a < Ami 



7TJ 

s — a i" 00 



ds' 



vr J Ami (s' — a)(s' — s) 



I s' — Am 2 t 1 \ s — a 
: , Ms) = 



n = 1,2. 



(20) 



(21) 



ds 1 



'Ami — s' 



7r j-00 (s' — a)(s' — s) V — s' 



f22) 



One can see from ()21|) that continuation of the loop to the second Riemann sheet is performed 



as: 



j(2) = -fW(s)-FW(o)- j L(s) = 
= JW - 2F« = 
= -J (1) -2(F(a)+L( S )). 



(23) 



The first expression seems the most convenient in numerical calculations. 
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